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We discuss the effect of mode transitions on the current (7) and temperature (T) dependent 
linewidth (A/) in nanocontact based spin torque oscillators (STOs). At constant 7, Af exhibits an 
anomalous temperature dependence near the mode transitions; A/ may either increase or decrease 
with T depending on the position w.r.t. the mode transition. We show that the behavior of A/ as a 
function of 7 can be fitted by the single mode analytical theory of STOs, even though there are two 
modes present near the mode transition, if the nonlinear amplification is determined directly from 
the experiment. Using a recently developed theory of two coupled modes, we show that the linewidth 
near mode transition can be described by an "effective" single-oscillator theory with an enhanced 
nonlinear amplification that carries additional temperature dependence, which thus qualitatively 
explain the experimental results. 

PACS numbers: 85.75.-d, 76.50.+g, 72.25.-b 



I. INTRODUCTION 

A spin-polarized current traversing a thin magnetic 
layer can exert a significant torque on the magnetiza- 
tion through the spin transfer torque (STT) effect 
The effect can be described as negative damping, lin- 
early proportional to the spin-polarized current, which 
at a certain threshold can overcome the natural Gilbert 
damping in the magnetic layer, allowing for coherent, 
large amplitude, excitation of spin waves. If the mag- 
netic layer is part of a structure with magnctorcsistance, 
such as a spin valve (SV) or a magnetic tunnel junction 
(MTJ), the excited spin waves can be used to generate a 
current- and field-tunable microwave voltage signal; the 
resulting device is commonly called a spin torque oscil- 
lator (STO). 8 Interest in STOs for microwave applica- 
tions is steadily increasing, due to their attractive com- 
bination of very large frequency tuning ranges^— effi- 
cient spin- wave emission in magnonic devices ^2r— very 
high modulation ratesji£~— sub- micron footprints^ and 
straightforward integration with semiconductor technol- 
ogy using the same processes as magnetoresistive random 
access memory i 25 i 26 

A minimal spectral linewidth, A/, of the microwave 
signal is highly desirable for applications. While a num- 
ber of recent experimental studies have addressed the 
temperature dependence of Af in nanopillar STOs— ~— 
the study of the temperature dependent linewidth in 
nanocontact STOs is limited to a recent work by Schnei- 
der et. al£2- The theory of the origin of STO lincwidths 
and their temperature dependence is now well established 



for single spin-wave modest - — A key result is the strong 
impact that limited amplitude noise can have on the STO 
phase noise, via the strong amplitude-phase coupling. 
Gaussian (white) amplitude noise is transformed into col- 
ored phase noise, and the intrinsic Lorentzian line shape 
expected for an auto-oscillator with zero amplitude-phase 
coupling changes into a convolution of Lorentzian and 
Gaussian line shapes4£ The coupling also leads to a 
substantial enhancement, or amplification, of the ther- 
mal broadening, and can also lead to asymmetric line 
shapes near threshold^ The degree of coloring should 
also change with temperature, leading to a crossover from 
a linear temperature dependence of Af at low tempera- 
ture, to a square root dependence at high temperature.— 

All temperature dependent studies to date show tem- 
perature regions with unexpected behavior. In Ref. l29l . 
Af in the subthreshold regime narrows by a factor of 6, 
from 1.2 GHz to 200 MHz, when the temperature is raised 
from 20 K to 140 K. In Ref. [H, the slope of the temper- 
ature dependence even changes sign multiple times as a 
function of drive current, and is close to zero at the small- 
est Af. In Ref. H3, Af increases exponentially above a 
certain temperature; the concept of mode hopping was 
introduced to explain and model this dependence. The 
origin of these rather complex temperature dependencies 
is yet to be explained. More recently, a linear behavior 
of linewidth is observed for a certain range of temper- 
ature in magnetic tunnel junction based STOsj 31 ' 32 A 
saturation of linewidth is observed in both these studies 
for temperature below 100 K, which is not explained by 
the existing theories. In addition, the temperature de- 
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pcndcncc of power restoration rate observed in Ref |32| 
can not be explained by the single mode theory^ Thus 
the details of the temperature dependence of linewidth 
in STOs is far from being understood. 

In this work, we present a detailed study of the 
temperature-dependent linewidth in nanocontact STOs. 
While all measurements were carried out at current and 
magnetic field values where only propagating spin waves 
were generated) 12 ' 41 we found a large number of mode 
transitions as a function of current (at a fixed temper- 
ature T) and temperature (at a fixed current /). The 
measured linewidth is highly nonmonotonic both as a 
function of current and of temperature, with large en- 
hancements at currents or temperatures where mode 
transitions occurred. We show that the linewidth is 
very well fitted by the single oscillator theor} ' 37 ' 38 , if 
the so-called amplification factor is obtained directly 
from measurements. While this agreement is similar 
to that of Refs. [l9| and [29|, we find the temperature 
dependence of the linewidth does not agree with that 
obtained directly from calculations using the nonlinear 
single-oscillator theor y 37 ' 38 , from which typically a lin- 
ear dependence on T is obtained for the systems un- 
der study here. These observations indicate that the 
central mechanism for linewidth broadening in nonlinear 
single-oscillator theory applies here, too: The linewidth is 
driven by phase noise amplified by the coupling through 
the nonlinear frequency shift to power amplitude fluctu- 
ations. However, our results indicate that this coupling 
may itself have a nontrivial temperature (and current) 
dependence, especially near mode transitions. We will 
here show that extending the nonlinear single-oscillator 
theory to include two coupled modest leads to additional 
couplings between the phase and power fluctuations. Un- 
der some simplifying assumptions, these couplings lead 
to a changed power restoration rate and the final re- 
sult for the linewidth looks very much like that from 
the nonlinear single-oscillator theor y 37 ' 38 , but with an 
enhanced nonlinear amplification that carries additional 
temperature dependence. This explains qualitatively the 
observed temperature dependence of the linewidth near 
mode transitions. 



II. EXPERIMENT 

The results presented in this work are from a sin- 
gle nanocontact STO device with an e-beam pat- 
terned 50 x 150 nm 2 elliptical nanocontact fabricated 
on top of a 8x26 /xm 2 pseudo-spin- valve mesa based 
on Co 8 iFeig(20 nm)/Cu(6 nm)/Ni 80 Fe2o(4.5 nm), as de- 
scribed in Ref. 0. While not shown here, other nanocon- 
tacts of varying sizes were also studied as a function of 
temperature, and gave the same qualitative results. 

The experimental circuit is similar to that employed in 
Refs. fTol and l44l The signal generated from the STO was 
amplified using a broadband +22-dB microwave ampli- 
fier, and detected by a 20 Hz-46 GHz Rohde & Schwarz 
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Figure 1. (Color online) (a) Two-dimensional power spectral 
density map of / versus / at a magnetic field of jj.oH—1 T, 
applied at an angle of 80° to the film plane. Top inset 
shows two examples of mode transitions at 7=30.2 mA and 
35.3 mA respectively, where the left spectrum has two clearly 
resolved Lorentzian peaks, and the right spectrum shows a 
single broader, asymmetric peak that can still be well fitted 
by two Lorentzian functions. The bottom inset shows the in- 
verse power 1/p vs current and a linear fit (solid line), (b) 
Experimentally measured (red triangles) and calculated Af 
(blue solid line). The black dashed line represents a linear fit 
to linewidth using Eq.(l) for subthreshold currents. 



FSU46 spectrum analyzer. The measurement was per- 
formed in the default mode of the spectrum analyzer, 
mode for spectrum analysis, the so-called analyzer mode. 
We use a resolution bandwidth of 10 MHz and video 
bandwidth of 10 kHz. The spectra were measured in 
the frequency range 13-25 GHz with a sweep time of 
100 ms. We also average 20 traces resulting in a total 
measurement time of about 6.4 s. The dc bias current 
is fed to the device by a current source through a 0-26 
GHz bias tee connected in parallel with the transmission 
line. The temperature of the sample was varied in the 
range 300-400 K through use of a heating foil underneath 
the sample. Each measurement temperature was main- 
tained with a precision of 0.1 K using a thermocouple 
attached to the bottom of the sample and a software- 
based PID controller. All measurements were performed 
in a /io-ff=T T field applied at an angle of 80° w.r.t. the 
film plane. In this geometry only a propagating spin wave 
mod e 12,14 ' 45 ' 46 is excited, and the output power is close 
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Figure 2. (Color online) Map of frequency of the strongest 
mode versus temperature and bias I, showing the mode tran- 
sition with temperature. The solid lines are linear fits to the 
threshold current for the mode transitions Ii and 1% versus 
temperature. The dotted lines are the positions at which the 
behavior of the linewidth is discussed in Fig. [3] 
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III. RESULTS 



Figure 3. (Color online) Measured linewidth versus tempera- 
ture at (a) 28.4 mA, (b) 29.4 mA, and (c) 30.3 mA. The solid 
black circles (respectively the solid blue squares) denote the 
mode excited below (above) 72=30 mA at room temperature, 
(d) Integrated power versus temperature at 28.4 mA (solid 
black circles), 29.4 mA (open symbols), and 30.3 mA (solid 
blue squares). The dashed lines serve as visual aids. 



Figure [T] shows the current (I) dependence of the STO 
frequency at room temperature. In addition to the ex- 
pected linear blue shift with /, a large number of discon- 
tinuous jumps and other nonlinearitics can be observed. 
We argue that all these nonlinear features are related to 
mode transitions, some large, where two distinct peaks 
can be observed on the spectrum analyzer [the left spec- 
trum in the inset of Fig. HJa)], and others small, where 
only a single peak is observed, though with a significant 
increase in both nonlincarity and linewidth [the right 
spectrum in the inset of Fig. QJa)]. Similar mode tran- 
sitions have been observed in the literatur o 27 ' 47 ' 48 and 
numerical simulations have reproduced this behavior for 
in-plane fields42 

The mode transitions have a significant impact on Af 
vs. I, as shown in Fig. Q] (b). We define Af as the full 
width at half maximum (FWHM) obtained by fitting a 
single Lorcntzian function. In the case of two modes, we 
use the linewidth of the strongest mode (the mode with 
the highest output power). In the subthreshold regime, 
A/ decreases linearly with increasing /, which we at- 
tribute to the narrowing of the natural ferromagnetic res- 
onance (FMR) linewidth under the influence of the neg- 
ative damping associated with spin torquei 35 ' 38 At every 
mode transition position, we also observe a dramatic in- 
crease in A/ leading to a highly nonlinear dependence on 
/. It is noteworthy that a strong mode transition, and 
the associated increase in Af, can also be observed well 
inside the subthreshold regime, at about 25 mA. The ex- 
istence of mode transitions is hence not limited to states 
of steady precession, as in Ref. H^. 

In order to show the effect of temperature on mode 
transitions, we plot a map of measured frequency vs tem- 
perature and current, as shown in Fig. [5] At room tem- 
perature these transitions are located at about 1\=17 mA 



and -?2=30 mA. As T is increased, both ii and I2 move 
to lower values following a linear dependence (the solid 
lines in Fig. [2]). This T dependence of I\ and I2 has di- 
rect consequences for A/(T). To illustrate this, we have 
chosen three current values, shown by the dashed lines in 
Fig. [21 which lie below, on top of, and above the second 
mode transition. Figures [3] (a)- |3|c) show Af vs. T at 
these three currents, which clearly exhibit three dramati- 
cally different T dependencies: i) at 28.4 mA, we observe 
a nonlinear increase of Af with T, ii) at 29.4 mA, we ob- 
serve a nonmonotonic T dependence, and iii) at 30.3 mA 
we observe an nonlinear decrease in Af with T . It is 
quite obvious that none of the measured curves in Fig. [3] 
follow either a linear or a square-root T dependence, as 
expected from the theories of thermally induced phase 



Now we will compare our results with the single mode 
analytical thcoryi 35 i 38 According to this theory, Af of a 
nonlinear oscillator is given by 

A/ = r g (l-£), for/«/ th (1) 

=A/ L (l + l / 2 ),for I» J th , (2) 

where T g is the natural FMR linewidth, I the bias 
current, 7 t h the threshold current, and the nonlinear 

/ \ 2 

linewidth amplification is (1 + v 2 ) = 1 + f ^— J , where 
= ^ is the nonlinear frequency shift, and T p is the 
power restoration rate (r" 1 is the correlation time of the 
power fluctuations); A/l = r g E k ^ g ? is the intrinsic ther- 
mal linewidth, i.e., the linewidth of a linear (v = 0) os- 
cillator. Here, E(po) is the total energy of the oscillator. 
Above threshold (I ^> 7 t h), the nonlinear amplification 
of the linewidth is controlled by the ratio of the non- 
linear frequency shift N to the power restoration rate 



4 



r o . The reason for this 3 ^ is that power fluctuations 
couple to phase fluctuations through the nonlinear fre- 
quency shift N, and the linewidth is dominated by phase 
fluctuations. The linewidth increases when N is large, 
so that small power fluctuations give rise to large phase 
fluctuations, or if T p is small, so that power fluctuations 
remain for a long time during which they affect phase 
fluctuations. For the nanocontact under study, the non- 
linear damping Q is small 3 ^, and we can approximate 

{l + v 



1 



r„ di 



We first compare our experimental results for fixed T 
with theory 35 ' 38 In order to do so, we need to extract T g . 
We fit the initial decrease in linewidth with Eq. (p}, and 
obtain T g = (500 ± 20) MHz and I th = (29 ± 1) mA, as 
shown by the dashed line in Fig. Q] (b). Next, from the 
measured / vs /, we obtain df jdl and directly calculate 
the nonlinear amplification factor (1 + u 2 ), and find from 
a fit to Eq. © that A/ L - 67 kHz, for I > I th . This 
value of A/l corresponds to kT/E(p ) ~ 1.5 x 10~ 4 . As 
shown in Fig. Q] (b), the calculated A/ shows very good 
agreement with the experimentally measured linewidth, 
and also reproduces the dramatic increase in A/ which 
occurs around each mode transition. The agreement in- 
dicates that the nonlinear amplification of the linewidth 
is controlled by the nonlinear frequency shift N oc df / di , 
while the power restoration rate T p is constant. The 
agreement is lost for I <27 mA, as expected for currents 
below threshold ! 35 ' 38 We have also used the inverse power 
method 5 -^ to determine the threshold current as shown in 
the inset of Fig. [T] (a) . A fit of this data for current be- 
low 25 mA is shown by the solid line. From this fit it 
appears as if the STO is close to auto-oscillation already 
at about 24.5 mA, but gets interrupted by one or more 
subthreshold mode transitions. It is only at about 27-28 
mA that robust auto-oscillation begins. 

Next, we want to compare the temperature dependence 
(at fixed I) of A/l and (1 + v 2 ) as obtained from the ex- 
periment with theoretical predictions 3 - 5 -^. According to 
the theory, A/l(T) should be proportional to T, since 
it is the linewidth of a linear oscillator in contact with 
a thermal bath, while (1 + v 2 ) has a monotonic tem- 
perature dependence. Using the agreement between the 
calculated and measured lincwidths in Fig. [TJ we can now 
extract A/l and its temperature dependence, as shown 
in Fig. |U[a). Since the determination of (1 + v 2 ) is more 
accurate in regions between mode transitions, i.e., where 
df /di do not diverge, we use the average value of A/l 
for 30.5 mA< / < 31.5 mA, which excludes any mode 
transitions and is above threshold at all temperatures. 
A linear increase in A/l with T is observed. The solid 
and dashed lines are calculations based on the classical 
quasi-Hamiltonian formalism for spin waves ; 35 ' 36 ' 38 ' 51 ' 52 
which shows reasonable agreement with the experiment 
and also predicts a linear behavior similar to experiment 
even with the inclusion of the temperature dependence 
of Mo in the calculation (red solid line) . This calculation 
assumes single mode excitation but considers the nonuni- 
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Figure 4. Temperature dependence of (a) extracted linear 
contribution of linewidth A/l (solid symbols) and (b) the 
nonlinear amplification, (1 + v 2 ) (solid and open symbols). 
The solid red lines are calculation based with inclusion of 
temperature dependence of M a , where as the dashed blue lines 
are calculation assuming no temperature dependence of M B . 



form nature of propagating spin waves by " exchange nor- 
malization" of magnetic field, and normalization of vol- 
ume under the nanocontact ^ The parameters used are 
similar to those of Ref. [lU, the electron gyromagnetic 
factor: 7 = 1.76 x 10 7 rad/Oe, saturation magnetiza- 
tion: Mq(300K) = 640 emu/cm 3 , Gilbert damping pa- 
rameter: etc = 0.01, dimensionless spin-polarization effi- 
ciency: e = 0.2, the exchange length: A ox = 5 nm, and 
(I / hh)3Q0K = 5. The effective volume V e s is assumed to 
be 1.5 times that of the volume under the nanocontact. 
We use T g = (500 ± 20) MHz, as determined from the 
experiment. Calculation also predicts T g = 500 MHz for 
our experimental geometry. We note that the agreement 
with A/l with T was obtained only when I /Ith > 5. We 
attribute this to the fact that the analytical Eq. © is 
an asymptotic equation that is valid only for I » Ith ^ 
Basically we treated I /Ith as a fitting parameter, with 
the other parameters kept fixed at their reasonable val- 
ues, since the precise values of these parameters are a bit 
uncertain. 

In Fig.@Jb) we show the behavior of measured (1 + is 2 ) 
vs T (symbols) for the three current values of 28.4 mA, 
29.4 mA and 30.3 mA along with the calculated behavior 
for I /Ith = 5 (solid and dashed lines). The experimen- 
tal behavior of (1 + v 2 ) vs T is dramatically different for 
the three cases but very similar to the behavior of the 
linewidth as a function of T shown in Fig[3J In contrast, 
the calculations of single-mode theory predict a mono- 
tonic decrease of (1 + v 2 ) with T when the temperature 
dependence of M s is included. Hence the calculations 
agree with the experiment only for a limited range of 
temperature and when the STO is far from the mode 
transition. For example, at 28.4 mA (30.3 mA), (1 + v 2 ) 
is enhanced at higher (lower) temperature, which is close 
to the mode transition. Detail examination show that 
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Figure 5. Map of power (dB) vs frequency (/) and tem- 
perature T for the three example current values of 28.4 mA, 
29.4 mA, and 30.3 mA. The arrows indicate the presence of 
additional modes, the amplitude of which depends on tem- 
perature. 



this enhancement occurs when two modes arc observed. 



This can be clearly seen in Fig. [5] which shows the mea- 
sured power vs frequency (/) and temperature T for the 
three current values. These spectra show the presence 
of an additional mode (as shown by the arrows) for all 
there currents. The temperature dependence of the am- 
plitude of this additional mode has a clear correlation 
with the behavior of (1 + v 2 ) vs T in Fig. HJb). For 
example, the amplitude of second mode increases (de- 
creases) with temperature at 28.4 mA (30.3 mA). Thus 
our results indicate that when two modes are observed, 
the experimental (1 + v 2 ) is enhanced compared to the 
prediction of single mode calculation. 

IV. DISCUSSION 

We will now discuss the mechanism for the anoma- 
lous temperature dependence of the lincwidth. The ba- 
sic assumption is that in the presence of two mode, mode 
coupling near a transition can lead to an increase in the 
linewidth. The starting point is a set of coupled equa- 
tions for the complex amplitudes c.;, i = 1, 2 of the time- 
dependence of the modes,— 



I 

+ iuji (pi,p 2 )c 1 + [r + (pi,p 2 ) - T_ (pi,P2)\ ci - ke llp c 2 = 
^ + ilu 2 (pi,P2)c 2 + [r+ (pi,p 2 ) - T_ (pi,p 2 )} c 2 - ke llp ci = 0. (3) 



Here, T+ and T_ are the positive and negative damp- 
ing, and wi and w 2 the mode-frequencies; we have indi- 
cated the dependence of taj, r+, and T_ on the mode 
powers pi and p 2 . The equations contain a linear cou- 
pling term with complex amplitude ke lLp , with k real and 
k > 0. This term is not allowed on short time scales if 
uj\ 7^ io 2 . Here, however, we are interested in behavior 
over times much larger than the time-scale of the pe- 
riods of the modes or of thermal fluctuations. In that 
case, the coupling mediated through the linear coupling 
describes processes in which one mode can decay into 
the other through intermediate states and energy that is 
absorbed or released into other magnetic modes or a ther- 
mal reservoir. Such a process becomes more likely as the 
mode frequencies approach each other, with a concomi- 
tant increase in k. The experiments show a significant 
current and temperature dependence of the main mode 
frequency. Therefore, the linear mode coupling also has a 
strong current and temperature dependence, k = k(I,T) 
and (f = tp(I,T). In particular, the magnitude of k has 
maxima at currents and temperatures at which mode 
transitions occur. We will see that this coupling plays 
a key role. 

We now make some simplifying assumptions. First, we 
assume that the mode frequency uii only depends on pi 
and not on pj, j ^ i. Next, we assume that the system is 



close to, but above, threshold (recall that the threshold 
current is about 27 mA and the relevant current values 
are around 30 mA). We then expand Eq. ((3]) near Cj = 
and write the equations in terms of power amplitude and 
phase, Ci = % e -»(^i,o*-yt) where o is the threshold 

mode frequency. This leads to the following equations for 
the time dependence of the amplitudes Qi and phases (fi : 

^ = T g (I/I th - 1) Q 1 - (QQ 2 + PQl) Qi 



+kQ 2 J — - cos(<p - ip 2 + ipi) (4a) 

V ^2.0 

^ = T 9 (J/J* - 1) Q a - (QQj + PQl) Q 2 

+kQi> -^-cos(<p- ip 2 - ipi) (4b) 

V w i.o 

d -^ = -N lQ l + k9l^^ + ^-^) (4c) 
at Qi y cj 2 ,o 

^ = -N 2 Ql + k9±^M<P ~ ^ + W ). (4d) 
at Q 2 y u>ifi 

Here, is the nonlinear frequency shift, the thresh- 
old current, and Q and P the diagonal and off-diagonal 
nonlinear damping coefficients, respectively. We will for 



6 



simplicity assume that Ni = N2 = N . 



Next, we introduce the transformations^ Qi = 



V / P ( 



e+ir/2 



and Q 2 



r- ■ ( e+n/2 



whcrc p is the 



total power in the two modes. These transformations re- 
cast the description of the mode amplitudes Qi in terms 

I 



of p and 9, where p is total power and 9 describes how the 
power is distributed between the two modes. Inserting 
these in Eqs. (j^a)) and (|4b[) and assuming that the aver- 
age power pq is stationary and writing p = po + 5p, with 
dpo/dt = and Sp the power fluctuations, we obtain the 
following linearized equation for the power (linearized in 
power fluctuations about the average power po): 



d(p + Sp) 
dt 



2 (I/I th - 1) T g (p + Sp) - 2Q (pg + 2 Po S P ) - (P - Q) cos 2 9 (p 2 + 2 Po Sp) 



+k (p + Sp) cos t 



— ^ cos (<p + ip) + , — — cos (<p - ip) 

k>2,0 V w l,0 



(5) 



where ip = ip2 — <pi, with a time evolution given by 

dip 



at - a atx ■ a 1 1 1 - sin 61 / w 2,o . , M , l + sin(9 /cJi )0 . , , 

— = — iVpo sin £7 — JV opsin (/ + fc . / sml^ — ip) — k — . / sinfyj + ip) 

dt cos 9 V uji cos £? y cj 2 



(0) 



r 



Ignoring fluctuations for the moment, and keeping in 
mind that the power po is constant, the equations ([5]) 
and §Q describe a two-dimensional dynamically driven 
system in (#, V>)-spacc. The system under consideration 
here has, far away from a mode transition so that k w 0, 
a single stable fixed point 9 = —n/2 (9 = n/2) with 
all power in mode oji (002) well below (above) the mode 



transition. Near or at the mode transition, the system 
may have stable fixed points, unstable fixed points, or 
limit cycles. In either case we will assume that the ex- 
perimental linewidth arises from fluctuations in the total 
power and phase difference, and wc will therefore ignore 
fluctuations in 9. By enforcing the stationarity condition 
dpo /dt = we obtain from Eq. (O 



2 (I/I th - 1) r s - 2Q Po -(P-Q) po cos 2 9 + k Po cos 9 



^(cos(v? + ip)) + . /^£(cos(v? - ip)) 
^2,0 V w i,o 



0, (7) 



r 



where (. . .) denotes a suitable time-average over times 
long compared to the time scale of fluctuations (e.g., a 
limit cycle). Inserting this into Eq. ([5]), and separating ip 
into a regular part describing the slow time evolution 
of the phase difference of the two modes, and fluctua- 

I 



tions Sip, ip = ^ + Sip, and replacing cos-0 (sin?/') with 
(cos 'I') ( (sin fy) ) we obtain the following linearized equa- 
tions relating the fluctuations in power and phase angle 
difference: 



^ = 2 (I/I th - 1) T g Sp - AQ Po Sp + 2 (P-Q) cos 2 



-kpo cos 9 Sip 



^ {S in(v+n-J^{sH<P-n 

0J2,0 V w l,0 



kSp cos ( 



^(cos(^ + *)) 



^(cos(y>- *)) 

Wl,0 



(8) 



r 



and 



dSip 
~df 



-NSp sin 9 
1 — sin ( 



-kSip 
-kSip 



cos u 
1 + sill I 
cos 9 



^(cos((^ - #)) 

Win 



with "J satisfying 



dt 



^(cos(<^ + *)) (9) 

w 2,0 



-Npo sin ( 



+k 
-k 



sm ( 



cost^ 
I + sin ( 



cos( 



w 2,0 . / T N 

— — sm(ip — W ) 



^ si n(^ + *). 

^2,0 



(10) 
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We pause for a moment to note that Eqs. |5J to 
(fTU)l restricted to a single mode (fc = and cos 6* = 
0) are precisely the results of Kim, Slavin, and 
Tiberkevich2&2&£i£2, with T p = - (I/I th - 1) T g + 2Qp, 
describing the power fluctuations in the oscillator, and 
how the power fluctuations couple to the phase fluctu- 
ations through the nonlinear frequency shift N . It is of 
course this latter coupling that gives rise to the enhanced 
linewidth through the enhanced phase fluctuations. As 
we noted earlier, in the low-temperature limit, applica- 
ble here, the single-oscillator linewidth enhancement is 
described by the ratio of the nonlinear frequency shift 
N to the power restoration rate T p : Power fluctuations 
couple to the nonlinear frequency shift, and the longer 
the decay time of power fluctuations is {i.e., smaller r p ), 
the more power fluctuations can affect phase fluctuations. 
For the system under consideration here, Eqs. (JSJ and ([9]) 
show that the mode-coupling k leads to additional cou- 
pling between power and phase fluctuations. In general, 
the solutions to these equations, especially in the pres- 
ence of thermal fluctuations, are complicated. We can, 
however, gain some insight by assuming that ojip ps W2,o 
with W2.0 > w i,o an d consider the system far from a mode 
transition so that k is small and 9 = —tt/2 + S, with 
j < 1, and ip small and negative. For the nanocon- 
tact STOs, the nonlinear amplification N is large, and 
the nonlinear damping small. First, with N large and 
k small, we can neglect the terms in Sip on the right- 
hand side of Eq. ([9]) . This means that power amplitude 
fluctuations couple to phase fluctuations through N just 
as for the single oscillator. It follows that if the power 
amplitude fluctuations are enhanced or prolonged by the 
mode coupling so that 5p is enhanced or T p reduced by 
the mode coupling, Second, for TV large and the nonlin- 
ear damping small, at the fixed point 9 ps —tt/2 we have 
cos(<p — $) PS and sm(ip — ps —I. If we neglect the 
terms in Sip on the right-hand side of Eq. ((HJ, the net 
effect under these assumptions is to change the power 
restoration rate r p — > T p — kcos(9) sin(|</?|) with a con- 
comitant enhancement of the nonlinear amplification and 
the linewidth as the coupling term v between power am- 
plitude and phase fluctuations is given by v = Npo/T p in 
single mode theory^ This explains qualitatively why the 
observed dependence of the linewidth on temperature in 
general does not agree with the theoretical expression^ 
(Fig. 01). In the latter, the temperature dependence is 
driven by the stochastic thermal noise. In contrast, the 
experimentally determined nonlinear amplification con- 
tains a modified power restoration rate that includes the 
temperature (and current) dependence of k (and tp). 



V. CONCLUSIONS 

In conclusion, we have shown that the behavior of spin 
torque oscillator linewidths is to a large extent deter- 
mined by nonlinearities arising from a number of mode 
transitions. The mode transitions are observed at in- 
creasing current at fixed temperature, or at increasing 
temperature at fixed current. Near the mode transitions, 
the linewidth increases substantially. Nevertheless, both 
the current and temperature dependence of the linewidth 
are well described by the analytical single-oscillator the- 
ory using the nonlinear amplification extracted from ex- 
perimental data. In contrast, the temperature depen- 
dence of the linewidth near the mode transitions does 
not agree well with the single-oscillator analytical theory 
if the nonlinear amplification is calculated directly from 
the theory. The experimental data showed the presence 
of an additional mode where the nonlinear amplification 
is enhanced near the mode transitions. We have argued 
that a temperature-dependent mode coupling leads to re- 
duction of the power restorations rate, and therefore an 
enhancement of the nonlinear amplification and of the 
linewidth, and that this at least qualitatively explains 
the anomalous temperature dependence of the linewidth 
near the mode transitions. These results are important 
for the understanding of linewidth in spin torque oscilla- 
tors. 
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